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Abstract. Constant pressure Monte Carlo simulations of a coarse grained off- 
lattice model for monolayers of amphiphilic molecules at the air/water interface 
are presented. Our study focusses on phase transitions within a monolayer rather 
than on self aggregation. We thus model the molecules as stiff chains of Lennard- 
Jones spheres with one slightly larger repulsive end bead (head) grafted to a 
planar surface. Depending on the size of the head, the temperature and the 
pressure, we find a variety of phases, which differ in tilt order (including tilt 
direction), and in positional order. In particular, we observe a modulated phase 
with a striped superstructure. The modulation results from the competition 
between two length scales, the head size and the tail diameter. As this mechanism 
is fairly general, it may conceivably also be relevant in experimental monolayers. 
We argue that the superstructure would be very difficult to detect in a scattering 
experiment, which perhaps accounts for the fact that it has not been reported 
so far. Finally the effect of varying the chain length on the phase diagram is 
discussed. Except at high pressures and temperatures, the phase boundaries in 
systems with longer chains are shifted to higher temperatures. 



I. INTRODUCTION 

Amphiphilic molecules on water often form monomolecular layers (Langmuir monolay- 
ers): The nonpolar tails prevent them from dissolving in the water, and the water soluble 
head groups cause them to spread on the water instead of agglomerating into droplets. Such 
monolayers are of interest for various reasons: As model systems which can teach about the 
structure of lipid bilayers, omnipresent in every living organism; as model systems which 
display all kinds of ordering phenomena in two dimensions. Amphiphile monolayers on 
solid substrates (Langmuir-Blodgett films JT|) are studied because of their possible use in 
electronic and opto-electronic devices. 

If the tails are sufficiently long, Langmuir monolayers show a rich spectrum of phases 
IHH, depending on the temperature and the area density of the molecules. Among these 
are a disordered fluid state (liquid expanded) and numerous condensed phases with varying 
degree of intrinsic order. For example, the low-temperature phases are crystalline, and those 
at higher temperatures are hexatic ("liquid condensed"): The positional order is lost, but 
the orientations of the bonds connecting nearest neighbors have long-range correlations [|J . 
The chains in the low-pressure phases are often collectively tilted in one direction, whereas 
they stand on average upright in the high-pressure phases. Other types of order exist as 
well. The characteristics of the phase diagrams are in many respect similar for lipids, fatty 
acids, esters and alcohols. Experimental results suggest some general rules: Adding CH2 



units to the hydrocarbon chains has the same effect as reducing the temperature. Playing 
with the size of the head group (e.g, by varying the pH of the subphase ||, or by exchanging 
the head group HQ]) basically affects the tilting transitions. Moreover, there seems to be a 
relation between the effective size of the head groups and the direction of tilt J?]]. We will 
come back to that point shortly. 

These general features are common to monolayers of very different molecules, which share 
only a few distinctive properties, such as the existence of a polar head and of at least one 
flexible nonpolar tail. Hence studies of idealized models seem to be a promising approach 
to explore the connection between the monolayer properties and the molecular structure. 
Popular coarse-grained models for amphiphilic molecules are, for example, Lennard- Jones 
chains @||: Beads of two different types are connected together with springs to form the 
nonpolar "tail" units and the polar "head units" units. Bead-spring models have been em- 
ployed to study self-assembling micelles @-13] properties of bilayers [13,14] and monolayers 
|T5Hl9|. If one is interested in the internal structure rather than in self-assembly, the latter 
are most simply modelled by chains attached to one head unit, which is confined to a surface. 
A monolayer model of this kind has first been investigated in simulations by Haas, Hilfer 
and Binder [f0)|,|l6] and later by us |T7|-|T9|. A related model of entirely stiff chains has been 
studied by Opps and coworkers |20[ . 



Collective tilt can be induced in a way which is closest to nature by making the heads 
slightly larger than the tail diameter. Even if heads and tail monomers are identical, chains 
sometimes tilt because of intermolecular packing effects - tilted chains lock better into each 
other. However, this tilting mechanism depends strongly on the internal chain structure and 
thus on the particular choice of the model parameters. The relation between chain stretching, 
chain packing, and chain tilting are far from transparent [16]. If the tilt is induced by a 
mismatch of head and tail size, on the other hand, many properties of the tilted state can 
be understood from simple geometrical arguments. For example, the relation between the 
effective head group size, the surface pressure, and the direction of tilt mentioned earlier 
can be rationalized as follows: If the heads are just slightly larger than the tails, the chains 
can optimize the distance to four neighbor chains in the state tilted towards next nearest 
neighbors, but only to two neighbor chains in the state tilted towards nearest neighbors. 
Tilt towards next nearest neighbors is thus more favorable. When the head group size 
reaches a certain value, the chains turn to tilt towards nearest neighbors. This is because 
the hexagonal lattice of the chains gets distorted in the direction of tilt. In the next nearest 
neighbor tilted state, two distances remain unaffected by the distortion. In the nearest 
neighbor tilted state, however, all six distances to nearest neighbors increase as a result, and 
the larger head groups can settle in happily in the larger cage. The argument can be cast 
into a simple mathematical formula ]^l|]. It predicts that chains with sufficiently large heads 
tilt towards nearest neighbors at low surface pressures, and towards next nearest neighbors 
at higher pressures. This is exactly what is found experimentally [El. 



In a recent paper [plj, we have presented extensive Monte Carlo simulations of a bead- 



spring model with a 10 % larger head bead attached to a planar surface. We have calculated 
the phase diagram in pressure-temperature space and in the space of molecular area vs. 
temperature. The model exhibited a disordered fluid and a number of condensed phases, an 
untilted one (U), one with tilt towards next nearest neighbors (NNN), and at low pressures 
one with tilt towards nearest neighbors (NN). Hence the model reproduces some of the 
important states of Langmuir monolayers and can be considered as a useful minimal model 
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for the study of such systems. The phase diagram in pressure-temperature space is shown 
in Fig. 0. The goal of the present study was to extend that work in order to gain a better 
understanding of the connection between the chain structure and the phase behavior. To this 
end, we have varied the size of the head bead and the chain length. Surprisingly, we found 
that for heads which are more than 14 % larger than the tails, a new phase emerges with a 
modulated superstructure and an average tilt direction which is intermediate between nearest 
and next nearest neighbors. Such a superstructure has not been reported experimentally; 
however, we will argue that it would probably be difficult to identify with the available 
experimental methods. We will characterize the modulated phase and present a detailed 
phase diagram for a system with 20 % larger head beads. Then we will study the effect on 
the phase transition of varying the chain length and comment on finite-size effects. Systems 
with very large heads shall also be discussed very briefly. 



II. THE MODEL 

We study a system of short chains of length N, with one larger head bead confined to 
move within a planar surface (The xy plane, z — 0). Beads are not allowed to enter the half 
space z < 0. They interact with a truncated Lennard- Jones potential 

■ ( (cr/r) 12 — 2{<j Irf + vA for r < c a 
\~L. ! (r)={ VW ' J ~ , (1) 

I for r > c a 

where the offset v c = 2c~ 6 — c~ 12 is chosen such that Vu(r) is continuous at the cutoff 
r = c a. As long as the cutoff parameter c > 1, the interactions have an attractive part. 
Here, we use c = 2 for the tail bead interactions. The head beads are made purely repulsive 
by choosing c = 1 for the head-head and the head-tail interactions. Furthermore, the head 
diameter is slightly larger than the tail diameter, er# > a. The effective diameter for the 
head-tail interactions is (an + cr) /2. Beads in a chain are connected by springs of equilibrium 
length do with the spring potential 



V s (d) 



!f d 2 s In [l - (d - d ) 2 /ds 2 ] for \d - d \ < d 



(2) 



oo for \d — do\ > d$ 

(FENE potential). The chains are made stiff with a stiffness potential 

V A = k A -(l-cos6), (3) 

which acts on the angle 9 between subsequent springs and favors 6 = (straight chains). 

We have performed Monte Carlo simulations of n chains grafted on a parallelogram with 
periodic boundaries in the xy direction, under conditions of constant spreading pressure n 
in the xy plane. In order to avoid internal shear stress, both side lengths L x and L y and 
the angle a of the parallelogram were allowed to fluctuate during the simulations. The 
Monte Carlo moves include single monomer displacements, and rescaling of all coordinates 
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such that L x , L y , or a changes. They are accepted or rejected according to a Metropolis 



prescription with the effective Hamiltonian [22 



H = E + UA- nNT ki(A) , (4) 

where E is the internal energy, II the applied spreading pressure, and A = L x L y sin a the 
area of the simulation box. During the simulations, we have monitored the "internal pressure 
tensor" 

nS = ^<E^> + ^^. (5) 

and checked that it is diagonal and identical to T15 a /3 as it should. Here the sum i runs over 
all monomers, a, f3 over the x and y coordinate, Fi denotes the force acting on monomer i, 
and 5 a p is the unit matrix. 

In addition to the simulations, we have also performed a low temperature phonon expan- 
sion. The zero temperature ground state, which minimizes the enthalpy is characterized 
by the minimum enthalpy H and by the ground state configuration {{ff }, L x , L®, a }. 
Note that it is continuously degenerate with respect to the simultaneous translation of all 
monomers in the x or y direction (or, alternatively, the translation of the of the simulation 
box). Hence one has / = 3nN— n+1 nontrivial degrees of freedom left. At low temperatures, 
we approximate the enthalpy by the harmonic expansion 

Hharm = #0 + ^ U M U, (6) 

where u = {{(** - rf )/a}, (L x - L° x )/a, (L y - L° y )/a, a - a } 
is a dimensionless vector of (supposedly small) deviations from the ground state, and 

Mij = - — — 

OUiOUj 

the matrix of second derivatives. It has / + 2 Eigenvalues e s , two of which are zero due 
to the translational invariance mentioned above. For a ground state with uniform tilt, M 
can be diagonalized without too much computational effort, since already a simple Fourier 
transform in the xy plane partly diagonalizes it. In the case of the modulated phase, we 
preferred to keep all the coordinates explicitly. Having obtained the Eigenvalues, one can 
rewrite (|B|) as 

1 / 

Hharm = H + - ^ £ s(s 5 (?) 

1 8=1 

where the sum now runs only over the nonzero eigenvalues, and ( s are the corresponding 
normal coordinates. Using this expression, one can calculate the partition function 
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Here (3 = 1/ksT is the Boltzmann factor, and the prime indicates that the trivial contribu- 
tion of translating the simulation box in the xy plane has been omitted in the integral. One 
obtains the free energy 

1 f 

{3G = -\nZ = (3H + - £ In (3e s + const. (8) 

2 s =i 

The approximation is justified as long as the enthalpy measured from simulations agrees 
with the prediction of the harmonic approximation 

H(T) ^H + f -k B T. (9) 

In our simulations, this was the case for k B T jt < 0.1. Note that even for temperatures 
T — > 0, the harmonic approximation is not applicable in every respect. For example, it 
yields the wrong thermal expansion coefficient p3| . However, this is no serious problem 



since we we use it mainly to calculate (3G at some point in phase space (n , T ) here. 
Given this reference value, one can determine G at other temperatures and pressures from 
simulations by thermodynamic integration 

G(H,T) = G(H ,To) + k B T jf {dH'^ - dT ' J^fT^ ( 10 ) 

as long as the path T from (n ,T ) to (n, T) does not cross a first order phase transition. 

The comparison of the free energies in the different phases allowed to pinpoint the dis- 
continuous phase transitions between condensed phases, even though some of them were 
highly metastable over large portions of the phase space. If one transition point is known, 
lines of first order transitions can also be traced with the Clausius Clapeyron equation 

511 _ AH 

6T~TAA' ^~ l > 

where AH is the enthalpy difference between the coexisting phases, and AA the area dif- 
ference. Applying the Clausius Clapeyron equation is equivalent to performing a thermo- 
dynamic integration along a very special path, the transition line. We have located phase 
transitions with both methods, by use of eqn. flllD , and by thermodynamic integration with 
different paths. The agreement was generally very good. 



The parameters of the model are chosen as in Ref. flT5|J16| and in our previous work [19 
d = 0.7cr, d s = 0.2a, k s = lOOe, and k^ = 10e. Unless stated otherwise, our results refer to 
systems of n = 144 chains with length N = 7. These systems decorrelate on average within 
200-1000 Monte Carlo steps (MCS), where one Monte Carlo step consists of Nn = 1008 
attempted monomer displacements and one attempt to rescale L x , L y , and a. The systems 
were usually equilibrated over 70.000 MCS, and data were collected every 500st configuration 
over a period of 200.000 MCS or more. 



III. RESULTS 
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A. The modulated phase 



The phase diagram for <jh = 11, i-e., heads which are 10 % larger than the tail diameter, 



has been discussed in detail in Ref. [19]. It includes a disordered phase (LE) and a number of 
condensed phases with or without uniform tilt in one direction (LC-U, LC-NN and LC-NNN, 
see Fig. [I]). Such phases have been found numerous times in experimental monolayers @||. 
On increasing the head size, we observe at low pressures a novel phase with a modulated, 
striped superstructure. Figure ^] shows the projection of such a configuration onto the xy 
plane. The modulation turns out to be very stable: Once a system has assumed a modulated 
state, it stays modulated even if the temperature is raised almost to the melting temperature, 
and up to pressures of II « 90e/cx 2 are applied, where the equilibrium phase transition to 
a uniformly tilted phase is long passed. On the other hand, uniformly tilted states are 
metastable deep in the region of thermodynamic stability of the modulated phase. Hence 
hysteresis effects are strong. At chain length N = 7, the modulated state is observed for 
head sizes between 1.14 ~ <t# ~ 1.27. (If the heads are made even larger, the tails do not 
form closed tilted layers any more, but aggregate into small micelle-like structures.) We have 
also found it in systems of chains with length N = 6 and N = 8. The width of the stripes 
increases weakly with increasing pressure. A reliable determination of the equilibrium stripe 
width is however difficult because of the periodicity imposed by the boundary conditions. 
For the special case on = 1.2a and N = 8, we have varied the system size from 64 to 256 
chains, and found that the most favored stripe width within the region of thermodynamic 
stability of the modulated phase is approximately four. 

A closer inspection of Fig. |2| gives some insight in the origin of the modulation: The 
chains arrange themselves in rows, and the heads (crosses) form a hexagonal lattice which is 
slightly tilted with respect to the rows built by the chains. By means of such a construction, 
the distance between the rows can be kept smaller than the distance between the heads, 
which is favorable from the point of view of chain packing. The modulation thus stems from 
a competition between two incommensurable lengths, the head size and the chain diameter. 
It is penalized by lines of "defects" at the boundaries of the stripes, where the heads are 
localized right between two rows. The chains there have to bend strongly so that at least 
their last segments can join the rows. Note that the transition from a modulated state 
to a uniformly tilted state involves a rearrangement of either the head lattice or the chain 
rows, which probably explains the strong hysteresis effects. From the above discussion, one 
might suspect that chain flexibility is needed to stabilize the modulated phase. However, a 



very similar structure has been observed by Opps et al in systems of entirely stiff rods [20 
Nevertheless, we expect that the chain stiffness will significantly affect the region of stability 
in phase space of the modulated phase. 

Experimentally, superstructures like the modulation are often identified by satellite peaks 
in the structure factor, defined as 



nN 2 



5^exp(igr$) , (12) 



where the sum runs over all monomers in the system. We shall briefly discuss the form of 
S(q) in the different ordered phases Pl,|T9l|. The structure factor of an untilted state has peaks 
in the xy plane at q z = which correspond to the reciprocal lattice of the hexagonal lattice. 



These Bragg rods |24| are very sharp in the x and y direction and broad in the z direction, 
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with a width which is inversely proportional to the layer thickness. In the uniformly tilted 
state, the Bragg rods move away from q z = on a plane which is perpendicular to the long 
axis of the chains 0. In addition, the positions of the rods in the x and y direction shift 
due to the distortion of the lattice. In the modulated phase, satellite peaks emerge in the 
direction perpendicular to the stripes. This is demonstrated in Fig. 

In grazing incidence x-ray studies of Langmuir monolayers, the situation is complicated 
by the fact that the diffraction patterns average over all domain orientations in the monolayer 
plane (powder average). Hence one only obtains information on the radial average S(q\\, q z ), 
with gy = yjq% + q^- Moreover, the intensity of the Bragg rods decays rapidly, so that only 
the lowest order peaks can be measured. The average structure factor S(q\\, q z ) in this g-range 
is displayed in Fig. § (c). For comparison, Figs. § (a) and (b) also show the corresponding 
pictures in phases with uniform tilt towards nearest and next nearest neighbors. In an 
untilted phase, the six first order Bragg peaks of the hexagonal lattice all collapse on one 
peak centered around q z = 0. In the tilted phases, they split up: Their q\\ positions differ 
due to the distortion of the lattice, and unless they belong to m vectors perpendicular to the 
direction of tilt, they move out in the q z direction. Note that for symmetry reasons, some 
of the peaks still collapse in the states with uniform tilt towards nearest and next nearest 
neighbors. In the modulated state, they are all separated from each other. However, this 
would equally be the case in a state with uniform tilt towards an intermediate direction 
between nearest and next nearest neighbors. Moreover, two of the peaks are so close that 
they can easily appear as just one peak if the experimental resolution is not extremely good. 
The satellite peaks practically vanish at these small gy-vectors. Hence the diffraction pattern 
of a modulated phase would most likely be interpreted as one of a phase with uniform tilt 
towards nearest neighbors, or as one of a phase with intermediate tilt direction. Such phases 
have indeed been reported. 

We conclude that if there is no experimental evidence for a modulated phase, there is 
also nothing that would exclude such a possibility. The modulation would be very hard to 
identify with a diffraction experiment. Considering how stable it is in our model over a wide 
parameter region, and considering that the underlying mechanism is a very general one - 
the competition between two incommensurable lengths characterizing head group packing 
and chain packing - we believe that there is a good probability that it might also be present 
in some real Langmuir monolayers. 



B. Phase behavior for chains with head size ah = 1.2 

We have seen that increasing the head size has at least one important effect on the 
phase behavior, it brings the modulated phase into existence. In order to study further 
consequences in detail, we have calculated the phase diagram for the case of chain length 
N = 7 and head size an = 1.2. Due to the large hysteresis effects mentioned earlier, 
the phase transition between the modulated phase and the uniformly tilted state had to be 
determined by thermodynamic integration methods at temperatures T < 0.8e//ce. At higher 
temperatures, hysteresis effects disappeared, but the difference between modulation and 
uniform tilt disappeared too, and the transition point became difficult to localize. For the 
other phase transitions between tilted and untilted, ordered and disordered states, hysteresis 
effects were negligible. These phase transitions leave their signature in various quantities: 
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Fig. |5] shows the area per molecule as a function of temperature for various pressures. 
At low pressures, one identifies one single phase transition accompanied by a large jump 
of the molecular area, whereas at high pressures, the curves indicate the presence of two 
singularities. The situation thus seems similar to that at head size ajj = 11 (Fig. Hj), 
where one observes one transition from a tilted ordered to an untilted disordered state at 
low pressures, and two decoupled transitions (one tilting transition and one order /disorder 
transition) at high pressures. 

To check that this interpretation is correct, we have determined the length of the average 
projection of the head-to-end vector of the chains on the xy plane 

R*» = (13) 

which is a suitable order parameter for collective tilt. Here [x] and [y] denote the x and 
y component of the head-to-end vector, averaged over the chains of a configuration, and 
(•) denotes the thermal average over all configurations. It is shown for various pressures 
in Fig. |]. As expected, it drops to zero at the temperature which corresponds to the 
first singularities in the area-temperature isobars. One can also inspect average tilt angle 
9 between the head-to-end vector and the surface normal (Fig. ^). Since this quantity 
does not measure azimuthal symmetry breaking, it is nonzero even in the "untilted" state. 
However, it exhibits a sharp kink at the tilting transition, which is in fact easier to localize 
than the kink in R xy . (Note that R xy also never reaches strictly zero due to the finite size 
of the system.) 

The quantity which measures the order /disorder transition is the hexagonal "bond ori- 
entational order parameter" ty§ of two dimensional melting, 



i n 6 

«r E E eMMjk, 

j=l k=l 



(14) 



The first sum j runs over all heads of the systems, the second k over the six nearest neigh- 
bors of j, and (pjk is the angle between the vector connecting the two heads and an arbitrary 
reference axis. As illustrated in Fig. |8|, the parameter drops to zero at the tempera- 
ture corresponding to the second singularity in 9 and A/n, indicating the transition to a 
disordered phase where nearest neighbor directions have no long-range correlations. 

The resulting phase diagrams are shown in the pressure-temperature plane and in the 
area-temperature plane in Figs. § and [K| At low temperatures, the modulated phase su- 



persedes the state with uniform tilt towards nearest neighbors. However, the transition 
pressure to the state with uniform tilt towards next nearest neighbor decreases with increas- 
ing temperature, whereas that of the (metastable) transition between two uniformly tilted 
states increases. Hence a region is spared out at higher temperatures where uniform tilt 
towards nearest neighbors is stable. Comparing the phase diagram Fig. | with that of Fig. 
|l|, one notices one more qualitative difference: The slope of the transition from tilted to 
untilted state is negative for the smaller heads o# = 1.1, whereas it is partly positive in 
the system with larger heads an = 1-2. On increasing the pressure at constant temperature 
T « 1.5e/ks, one observes reentrant behavior: First a transition from an untilted to a tilted 
state takes place, then a second one back to a tilted state. At first sight, it seems strange 



that the chains should start to tilt on being pushed together. Fig. [10] suggests a possible 



explanation: Tilt is obviously most stable, i.e., persists up to the highest temperatures, at 



S 



the area per molecule A/n = la 2 . If the chains are squeezed together to lower molecular 
areas, they stand up. On the other hand, if the chains are allowed to occupy more space, 
they are in less close contact with each other and disorder more easily. Thus they remain 
tilted individually (cf. Fig. 0), but they lose the collective tilt. 

The phase diagram was determined from Monte Carlo simulations of systems of 144 
chains. In order to get a rough estimate of the finite-size effects, we have performed a 
few simulations of systems of 900 chains at pressures II = le/er 2 and II = 50e/a 2 . Fig. 



IT] suggests that finite-size effects tend to shift the apparent melting transition to slightly 
higher temperatures in small systems. Yet the shift is so small at pressure II = 50e/cr 2 that 
the transition temperatures in the two systems can still be considered the same within the 
error. At the lower pressure II = le/a 2 , one has to keep in mind the additional complication 
that the ordered state is modulated. One would thus expect that some particle numbers are 
more "commensurate" with this superstructure than others, and systems of different size 
are not directly comparable. Moreover, the relaxation times are very long, and the large 
systems could not optimize the stripe width within the simulation time (400.000 MCS per 
data point). Hence it is perhaps not surprising that they were found to melt earlier than 
the smaller system. No systematic finite-size effects were found for the decoupled tilting 
transition at pressure II = 50e/a 2 (Fig. |i~2j). 

Finite-size effects do of course affect the apparent phase boundaries, and the actual 
positions of the boundary lines in Figs. |9| and |10] have to be regarded with caution. A deter- 
mination of the phase diagram using systematic finite-size scaling methods would obviously 
be desirable. However, this would not only involve the study of much larger systems, but 
also require a much better statistics of the data. The computational costs would increase 
by orders of magnitude. We believe that our simulations of a relatively small system and 
simulation runs of length 200.000 - 400.000 MCS are sufficient to reliably characterize the 
phase behavior in many respects. Of course, important questions remain open. For exam- 
ple, we are not able to distinguish between a solid state and a hexatic state in the ordered 
region. We have not fully characterized the transition from the tilted to the untilted state 
- whether it is first order or continuous, whether it is one transition or whether one has 
an intermediate "unlocked tilted" state, where the tilt direction is not coupled to the local 
chain lattice etc. These questions will have to be addressed in future studies. 



C. Effect of chain length 

Experimentally, it is often observed that increasing the length of the hydrophobic chains 
causes the phase transitions to shift to higher temperatures || . One possible explanation is 
that the internal energy per molecule scales roughly linearly with the chain length, and the 
entropy gain associated with a phase transition is approximately independent of the chain 
length. Whereas the first assumption seems reasonable, the second one is more questionable. 
It is conceivable for the case of phase transitions between condensed phases, which affect the 
orientational degrees of freedom of whole chains. However, the melting transition involves a 
considerable gain of conformational entropy within the chains, which is obviously not chain 
length independent. 

In order to assess the influence of the chain length on the phase behavior in our model, 
we have studied systems of chains with length N = Q,7 and N = 8 at the spreading pressure 
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II = le/cr 2 and II = 50e/cr 2 . The results for the average tilt angle 9 and the order parameter 
of melting are shown in Figs. 13 and 14. The transition temperature of the low-pressure 
coupled phase transition between the tilted ordered and the untilted disordered state at 
indeed shifts to higher temperatures with increasing chain length (Fig. [H|). The same holds 
for the decoupled tilting transition at higher pressure (Fig. [TJ] (a)). The melting transition 
between the untilted ordered phase and the disordered phase, on the other hand, moves to 
lower temperatures as the chain length is increased (Fig. [14] (b)). 

It is instructive to examine the way how the area per molecule depends on the chain 
length, Fig. [15]. At low pressure, the area per molecule decreases with increasing chain length 
at all temperatures - presumably an effect of the enhanced attractive chain interactions (Fig. 
|T5| (a)). This remains true at II = 50e/cx 2 in the tilted state. However, the situation turns 
around at the tilting transition, and the largest chains occupy the largest area at all higher 
temperatures. Thus the area per molecule at this high pressure is determined by the entropic 
repulsion between the chains rather than by the energetic attraction between the beads (Fig. 

mm 



IV. SUMMARY 

We have examined the phase behavior of a system of stiff short Lennard- Jones chains 
grafted on a two dimensional plane with variable head size. As a function of temperature and 
surface pressure, we found a disordered phase and a number of condensed ordered phases with 
or without collective chain tilt. The tilting transition and the melting transition are coupled 
at low surface pressure and decouple at higher pressures. In addition, we observe a new 
modulated phase for a range of head sizes, resulting from the competition between the head 
and the tail size. As we have argued, there is a possibility that such a phase is also present 
in real monolayers. It would be characterized by a tilt in a direction intermediate between 
nearest and next nearest neighbors (which is observed experimentally) and by satellite peaks 
in the structure factor at high q vectors (which are hard to resolve in scattering experiments) . 
With increasing chain length, most phase transitions are shifted to higher temperatures, 
except at high surface pressure the melting temperature. 

The model is versatile enough to provide a useful starting point for further investigations. 
Simulations of much larger systems would allow to study lattice defects and defect interac- 
tions, and to distinguish between a crystalline and a hexatic phase (i.e., one with long range 
bond orientational order, but only short range positional correlations 0]). This is of interest 
because Langmuir monolayers are among the few two-dimensional systems where there is 



strong experimental evidence for the existence of hexatic states |25 |. Since amphiphile layers 
often contain amphiphiles of different type, it would also be interesting to study mixtures of 
amphiphiles, e.g., long and short amphiphiles |26], or to add a few amphiphiles with chain 



'defects", i.e., built-in kinks in the chains |Z7|| . Other promising problems relate to the 
interactions between chains and large inclusions, e,g, model "proteins" p8p9f , and the way 
how these are affected by the vicinity of phase transitions. Furthermore, future studies will 
be devoted to a more refined treatment of the head groups. In particular, it will be necces- 
sary to somewhat relax the rigid constraint of confining them in an entirely flat surface - 
by allowing them to move also in the z direction and/or by allowing for undulations in the 
substrate. 



10 



ACKNOWLEDGMENTS 



We have benefitted from fruitful interactions with Harald Lange and Kurt Binder. Frank 
Martin Haas and Rudolf Hilfer let us have their simulation code, which has been the base for 
the development of the code used here. C.S. thanks the Deutsche Forschungsgemeinschaft 
(DFG) for support through the Graduiertenkolleg on supramolecular systems in Mainz, and 
F.S. thanks the DFG for a Heisenberg fellowship. 



[1] W. A. Barlow, Langmuir-Blodgett Films (Elsevier, 1980); M. C. Petty, Langmuir-Blodgett 
Films (Cambridge University Press, 1996). 

[2] V. M. Kaganer, H. Mohwald, P. Dutta, Rev. Mod. Phys., in press (1998). 

[3] G. M. Bell, L. L. Combs, L. J. Dunne, Chem. Rev. 81, 15 (1981); H. Mohwald, Ann. Rev. 
Phys. Chem. 41, 441 (1990); H. M. McConnell, ibid 42, 171 (1991); C. M. Knobler, R.C. 
Desai, ibid 43, 207 (1992). D. Andelman, F. Brochard, CM. Knobler, F. Rondelez, in Micelles, 
Membranes, Microemulsions and Monolayers, p. 559, W. M. Gelbart, D. Roux and A. Ben- 
Shaul eds. (Springer, 994). 

[4] B. I. Halperin, D. R. Nelson, Phys. Rev. Lett. 41, 121 (1978). 

[5] M.C. Shin, TM. Bohanon, J.M. Mikrut, P. Zschack, P. Dutta, J. Chem. Phys. 96, 1556 (1992). 

[6] M.C. Shin, TM. Bohanon, J.M. Mikrut, P. Zschack, P. Dutta, Phys. Rev. A 45, 5734 (1992); 
J. Chem. Phys. 97, 4485 (1992). 

[7] B. Fischer, E. Teer, C. M. Knobler, J. Chem. Phys. 103, 2365 (1995); E. Teer, C. M. Knobler, 
C. Lautz, S. Wurlitzer, J. Kildae, T. M. Fischer, J. Chem. Phys. 106, 1913 (1997). 

[8] B. Smit, A. G. Schlijper, L. A. M. Rupert, N. M. van Os, J. Phys. Chem. 94, 6933 (1990); B. 
Smit, P. A. J. Hilbers, K. Esselink, L. A. M. Rupert, N. M. van Os, A. G. Schlijper, Nature 
348, 624 (1990); J. Phys. Chem. 95, 6361 (1991); S. Karaborni, M. M. van Os, P. A. J. Hilbers, 
Langmuir 9, 1175 (1993). 

[9] B. Smit, P. A. J. Hilbers, K. Esselink, Intnl. Journ. Mod. Phys. C 4, 393 (1993); B. Smit, K. 
Esselink, P. A. J. Hilbers, N. M. van Os, L. A. M. Rupert, I. Szleifer. Langmuir 9, 9 (1993); K. 
Esselink, P. A. J. Hilbers, N. M. van Os, B. Smit, S. Karaborni, Coll. Surf. A 91, 155 (1994); 
S. Karaborni, K. Esselink, P. A. J. Hilbers, B. Smit, J. Karthauser, N. M. van Os, R. Zana, 
Science 266 254 (1994). 

[10] B. J. Palmer, J. Liu, Langmuir 12, 746 (1996); 12, 6015 (1996). 

[11] F. K. von Gottberg, K. A. Smith, T. A. Hatton, J. Chem. Phys. 106, 9850 (1997); 108, 2232 
(1998). 

[12] D. Viduna, A. Milchev, K. Binder, to appear in Macromol. Theory Simul., 1998. 
[13] R. Gotz, R. Lipowski, J. Chem. Phys. 108, 7397 (1998). 



11 



[14] D. Harries, A. Ben-Shaul, J. Chem. Phys. 106, 1609 (1997). 

[15] R. Hilfer, F. M. Haas, K. Binder, Nuovo Cimento 16, 1297 (1994); F.M. Haas, R. Hilfer, K. 
Binder, J. Chem. Phys. 102, 2960 (1995). F. M. Haas, R. Hilfer, K Binder, J. Phys. Chem. 
100, 15290 (1996). 

[16] F. M. Haas, R. Hilfer, J. Chem. Phys. 105, 3859 (1996). 

[17] F. Schmid, C. Stadler, H. Lange, Computer Simulations in Condensed Matter Vol. 10, pp. 37, 
D. P. Landau, K. K. Mon, B. Schiittler eds. (Springer, 1997); Coll. and Surf. A, in press. 

[18] C. Stadler, Dissertation Universitat Mainz (1998). 

[19] C. Stadler, H. Lange, F. Schmid, submitted to Phys. Rev. E (1998). 

[20] S. Opps, B. Yang, C. Gray, D. Sullivan, to be published; S. Opps, private communication 1998. 
[21] F. Schmid, H. Lange, J. Chem. Phys. 106, 3757 (1997). 

[22] M. P. Allen, D. J. Tildesley, Computer Simulation of Liquids, (Oxford University Press, 1987). 

[23] N. W. Ashcroft, N. D. Mermin, Solid state physics (Saunders College, 1976). 

[24] Strictly speaking, one has no "Bragg" rods in a strictly two-dimensional system like ours, since 
there is no true long-range positional order. However, for all practical purposes, the peaks are 
indistinguishible from delta-singularities. 

[25] C.A. Helm, H. Mohwald, K. Kjaer, J. Als-Nielsen, Biophys. Journ. 52, 381 (1987); K. Kjaer, 
J. Als-Nielsen, C.A. Helm, L.A. Laxhuber, H. Mohwald, Phys. Rev. Lett. 58, 2224 (1987). 

[26] D. Chowdhury, P. K. Maiti, S. Sabhapandit, P. Taneja, Phys. Rev. E 56, 667 (1997). 

[27] Y. K. Levine, A. Kolinski, J. Skolnick, J. Chem. Phys. 98, 7581 (1993). 

[28] A. Baumgartner, J. Chem. Phys. 103, 10669 (1995); Bioph. J. 71, 1248 (1996). 

[29] T. Sintes, A. Baumgartner, Bioph. J. 73, 2251 (1997); Physica A 249, 571 (1998). 



12 




12 3 4 

Temperature T 



FIG. 1. 

Phase diagram of a model of endgrafted bead-spring chains with 10 % larger heads in the 
pressure-temperature plane. The model is described in section ||. The spreading pressure 
IT is given in units of e/cr 2 , and the temperature T in units of e/ks- LE denotes disordered 
phase, LC-NN ordered phase with tilt towards nearest neighbors, LC-NNN ordered phase 
with tilt towards next nearest neighbors, and LC-U untilted ordered phase. The transition 
between LC-NN and LC-NNN could not be located at pressures above II = 20e/<r 2 . See text 
for more explanation. From Ref. jnj. 
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FIG. 2. 

Projection on the xy plane of a configuration in a modulated state. Crosses mark the 
positions of head beads. Parameters are II = le/a 2 , T = 0.1e/k B , o~h = 1.2a. 
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FIG. 3. 

Integrated structure factor S int ,(q\\) = JqO^^j^dQz S(q) in the modulated phase. For the 
clarity of the presentation, we have integrated over the Bragg rods in the (tilted) Bragg plane 
q\ = —q\\et/ cos0, where 9 is the tilt angle, the unit vector e t points towards the average 
tilt direction, and q\\ = (q x ,q y ,fy- The interval of integration was Aq z = 2ir/d cos9. One 
clearly discernes rows of satellite peaks between the main peaks, indicating the presence of 
a periodic superstructure. 
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FIG. 4. 

Powder averaged structure factor S(q»,q z ) (a) in a state with uniform tilt towards nearest 
neighbors (T = 0.1e//cs,n = le/a 2 ,an = l.lcr) (b) in a state with uniform tilt towards 
next nearest neighbors (T = 0.1e//cB,IT = 50e/cr 2 ,<7f/ = l.lcr) (c) in a modulated state 
{T = 0.1e/k B ,U = le/a 2 ,a H = 1.2a). 
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FIG. 5. 

Area per molecule A/n in units of a 2 vs. temperature T in units of e/ks at head size 
an = ^-2cr for different pressures II in units of e/a 2 as indicated. 
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FIG. 6. 

Order parameter R xy in units of a 2 vs. temperature T in units of ej ks at head size oh = l-2cx 
for different pressures II in units of e/a 2 as indicated. 
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FIG. 7. 

Average tilt angle 6 in degrees vs. temperature T in units of e/ks at head size an = l-2cx 
for different pressures II in units of e/a 2 as indicated. 
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FIG. 8. 

Order parameter \l/ 6 vs. temperature T in units of e/k B at head size a H = 1.2a for different 
pressures II in units of e/cr 2 as indicated. 
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FIG. 9. 

Phase diagram of a monolayer of chains with heads of size <jh = 1.2a and length N = 7 
beads. Units of pressure are e/a 2 , and those of temperature e/ks- LE denotes the disordered 
(expanded) phase, LC-U an untilted ordered phase, LC-NN and LC-NNN ordered phases 
with uniform tilt towards nearest and next nearest neighbors, respectively, and LC-NN mod. 
the modulated phase. The dotted line indicates the metastable phase transition between 
two uniformly tilted phases in the region of stability of the modulated phase. The transition 
between LC-NN and LC-NNN at high temperatures is washed out and hard to localize. 
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FIG. 10. 

Same as Fig. || in the area-temperature plane. Area per molecule A/n is given in units of 
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FIG. 11. 

Melting order parameter \l/ 6 for systems of 144 chains (solid line) and 900 chains (dotted line) 
as a function of temperature in units of e/ks- Pressures are II = le/a 2 (a) and II = 50e/a 2 
(b). Head size is a H = 1.2a and chain length N = 7. See text for discussion. 
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FIG. 12. 

Average tilt angle 6 for systems of 144 chains (solid line) and 900 chains (dotted line) vs. 
temperature in units of e/ks- Pressure is II = 50e/a 2 and head size on = 1.2a. 
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FIG. 13. 

Average tilt angle 9 (a) and order parameter of melting \l/ 6 (b) vs. temperature (in units of 
e/ks) at pressure II = le/a 2 for chains of length N = 6 (solid lines), N = 7 (dashed lines), 
and N = 8 (dotted line). Head size is <jh = 1-2(7. 
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FIG. 14. 

Average tilt angle (a) and order parameter of melting \l/ 6 (b) vs. temperature (in units of 
e/ks) at pressure II = 50e/a 2 for chains of length N = 6 (solid lines), N = 7 (dashed lines), 
and N = 8 (dotted line). Head size is an = 1.2o\ 
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FIG. 15. 

Area per molecule A/n in units of a vs. temperature in units of e/ks at pressure II = le/a 2 
(a) and IT = 50e/<r 2 (b) for chains of length N = 6 (solid lines), N = 7 (dashed lines), and 
N = 8 (dotted line). Head size is a H = 1.2a. 
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